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This note is based upon the following excellent textbooks: 
Jackson: Classical Electrodynamics 

Panofsky and Phillips: Classical Electricity and Magnetism 
Griffiths: Introduction to Electrodynamics 
Raychaudhuri : The Theory of Electricity and Magnetism 
Feynman Lectures , vol . 2 

You are always advised to read the original textbooks. Remember that the supplementary problems 
form an integral part of the course. 

I will use the rationalised Lorentz-Heaviside system throughout. That system is explained in 
Section 1. I will also assume, as prerequisite, that you know (i) the ordinary noncovariant form of 
Maxwell's equations; (ii) the basic postulates of the Special Theory of Relativity, and (iii) how to 
get the Lagrangian and the Hamiltonian of a particle. In fact, for a certain section, I have to use 
the Lagrangian formulation for a field, but probably you have already encountered that in your 
Classical Mechanics course. Even if you have not, the formulation is analogous. 

I will always consider electromagnetic fields in vacuum, i.e., E = D and B = H. In fact, I will 
never use D and H and associated quantities like the electric polarisation P or magnetisation M. 

A four-dimensional vector will be labeled by a Greek index; all Greek indices will run 
from to 3, being the time component. I will use the flat space-time Minkowski metric 
ilfj-u = diag(l, —1, —1, —1). All repeated indices arc implicitly summed over. 

This note is brief, and please remember that the problems are an integral part of the course; you 
must do them before proceeding to the next section. I could not discuss a number of interesting 
topics just for want of time. For example, I would like to discuss the physics of accelerated charges 
in more detail, including the frequency dependence of radiated power; also the theory of half- 
advanced and half-retarded potentials. There may be numerous typos. Please feel free to inform 
me about them. You can contact me either at my office, or at akphy@caluniv . ac . in, where I am 
always available. 
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1 Maxwell's Equations 



Why should we take this course? Is something exciting happening in classical electrodynamics? 
The answer, unfortunately, is no; it is an age-old and established subject, like classical mechanics, 
and almost all the fundamental discoveries were made in the 19th century. 

Like classical mechanics, it is a tool (well, it is classical mechanics, only the force is provided 
by electric or magnetic fields, and you study the motion of charged particles in these fields, or the 
properties of the fields themselves), which is applied in many branches of physics, most notably in 
astrophysics, plasma physics (either man-made or astrophysical) , and particle physics, and in a lot of 
applied branches as well. Let me give you just an example. You know that Special Relativity comes 
into play when the velocities are very large, comparable to that of light, and General Relativity is 
relevant for dealing with gravity. In astrophysics, you encounter both these situations. There arc 
all sorts of radiation in the sky, and a proper knowledge of them can help us unravel the mystery of 
the origin of the universe, or at least the galaxies. You can map the sky not only with visible light 
but also with invisible electromagnetic waves, from microwaves to gamma-rays. So you must know 
how and where these radiations are generated. Relativistic classical electrodynamics helps you to 
do that. 

But electrodynamics has been quantised, and quantum electrodynamics is known to be one of the 
greatest intellectual achievements of mankind. No other theory, not even gravitation (Newtonian 
or Einsteinian) , has been tested to such precision. Why, then, we still use the classical theory? 
Indeed, there are places {e.g., if you wish to calculate the scattering cross-section of an electron in 
the field of another electron, or the well-known Compton scattering) where a quantum calculation 
is easier and more precise. But these systems are microscopic; when you apply your theory to a 
macroscopic system (like a star), you have to use some kind of averaging over the ensemble, and 
come back to the classical regime, thanks to the correspondence principle. Also, the classical theory 
deals with concepts like electric and magnetic fields, which are easier to measure. 

A good point to start this course is the set of equations that tells you almost everything about 
classical electrodynamics (the force law of Lorentz completes the set). By the word classical we 
mean that there are no photons, only the electromagnetic wave. Electrons are the objects with 
which this classical field interacts. Sometimes the electrons are treated in a quantum-mechanical 
way; that is called a semiclassical approximation. Of course, we have no time here to go into a full 
quantum-mechanical theory of the interaction of electrons with photons. 

In most textbooks {e.g., Panofsky and Phillips) you will see the use of the SI, sometimes called 
the MKS, system. In that system Maxwell's equations (in vacuum) read 

V.E = ^, 

VxB = m + fio£Q-g^, 

V.B = 0. (1) 

The electric charge is expressed in coulomb; the charge of an electron is —1.6 x 10^^*^ coulomb. 

In Raychaudhuri, as well as in the second edition of Jackson, you will encounter the Gaussian 
system, which, I must say, is gradually going out of use ^: 

V.E = 47rp, 

^Even though the system has become almost archaic, I will still recommend Raychaudhuri's book to those who 
are not lucky to have him as a teacher — he died in 2005 — but would like to have a glimpse of the physics insight 
of a teacher sans pareil. 
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V.B = 0. (2) 

The Lorentz force law in these two systems is 

F = g(E + V X B) (SI), F = g + X B^ (Gaussian). (3) 

In this note we will use the so-called rationalised Lorentz- Heaviside (RLH) system. This system 
is characterised by a scaling of the fields and charge and current densities of the SI system: 

E57 = —;=^RLH, Bsi = -s/JIoBrlh, qSI = V^OQRLH, jsi = —=jRLH- (4) 

This seems to be a complicated scaling, and to top it all, we set c = 1. That is not a blasphemy. 
It just tells you that the units of length and time are related, and when you say the length of 1 
second you actually mean the length light travels in 1 second (so the length of one year is actually 
one light-year) . It also tells you that any velocity, apart from being a dimensionless quantity, must 
be a number whose magnitude is less than or equal to one. Some quantities get related too: energy 
and mass now have same dimension, so we can safely talk about a proton having a mass of 938 
MeV 2. 

But what do we gain? You may not believe it, but we have abolished all factors of /^o, eo; and 
c from our subsequent discussions. E and B now stand on the same footing, so do p and j. Of 
course, physics does not change an iota, and if you like, at the end of your calculation, you can go 
back to the familiar SI result with the help of eq. (4). And after rescaling, it is better to make a 
dimensional analysis; you may have to introduce suitable powers of c in the expression 

To verify this claim of elegance, let us look at the Maxwell's equations, with eq. (4) and the 
relation iiqEq = 1/c^: 

V.E = p, 

V, „ SB 

V.B = 0. (5) 

An important point to note is that if p = j = 0, the equations are symmetric under the interchange 
E ^ B, B ^ — E. We will come back to this later. 

Under the parity transformation x — x, it is clear that E — E, B ^ B, j ^ — j. Under time 
reversal t —t, E does not change sign — it cannot if it is created by a static charge configuration 
— but B — B, which is intuitively easy to understand: a magnetic field is created by a current 
and under time reversal, the current reverses its direction, so does B. 

Of course, there is no free lunch, but the price to pay is small. We cannot express the electric 
charge in the conventional unit of coulomb. Rather, the charge is something like 0.3 unit (this 
comes from the fine structure constant, something you will learn later, whose experimental value is 



^This is the system that is used by astrophysicists and particle physicists. They also use h = 1 which relates mass 
with length and time: [M] = [L]~^ = [T]~^. We won't need it for a discussion of classical electrodynamics. You may 
go further and use fcs — 1; that will relate temperature with mass. 

^Our system only abolishes no, eo and c, but not numerical factors like 47r, so there should not be any problem 
with dimensional analysis. An example will be given later. 
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approximately 1/137, and the theoretical expression, in the RLH system, and with ^ = 1, is e^/47r 
— in the SI system it is / At:£qTic). There are people who believe that this number is more handy 
than something like 10~^^! 

Q. B is derivable from a vector potential A: B = V x A. Show that for a uniform magnetic field 
B, one can write A = ^B x r. 

Q. Check eq. (5). How does the Lorentz force law look like? 

Q. With c = 1, show that the Lorentz transformation equations in 1 + 1 (i.e., one space and one 
time) dimension can be written as t' = tcoshO — xsmh.6, x' = xcosh^ — isinh^. How is 9 related 

with v7 

Q. With c = h = l, find how MeV is related to s~^ and fm~^ (1 fm= 10~^^ m). 
Q. Show that the energy density of the electromagnetic field and the Poynting vector are respec- 
tively given by i (E^ + B^) and E X B in the RLH system. 



2 Four-vectors 

There are two postulates of the Special Theory of Relativity (STR): physical laws are invariant 
in all frames which are mutually inertial, and the velocity of light in vacuum is a constant in all 
inertial frames. A consequence of the second postulate is that nothing can travel faster than light 
(in vacuum). This consequence rules out action-at-a-distance; every information must take some 
time to proceed from one point to another ^. So if I shake an electron here, another electron at 
a distance of one light-year will feel it not before one year, and if the sun vanishes right now, the 
earth will fly off in a tangent only after eight minutes (assuming that the gravitational information 
travels at the speed of light). The objects which carry these informations are called fields. Wc may 
even quantise these fields and get the corresponding field excitations, the particles. It is just a small 
matter that electrodynamics was quantised long ago (and some other forces which are analogous to 
electrodynamics were born quantised), and gravitation resisted all attempts (and by the smartest 
minds on this planet!) of quantisation. 

Now to the more mundane subject of four-vectors. A three-vector has three spatial components, 
they satisfy certain transformation laws, and their products are defined in a certain way. Four- 
vectors have one temporal and three spatial components. They are written as 

A^' = {A^,A). (6) 

Conventionally, the Greek indices run from to 3 (and the Latin indices from 1 to 3, i.e., over the 
spatial components only). The zero-th component is the time component, and components 1 to 3 
are the usual spatial components. 

Af^ will be a four-vector if and only if its components transform like the transformation of x^, 
the position four- vector, defined as x'^ = (t,x) (remember c = 1): 

x'^ = A^x" (7) 

where is the Lorentz transformation matrix. For a boost along the x-direction, it looks like 
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^That's the problem with Coulomb's law. It just tells you the force between two charges, and is of an action-at- 
a-distance form. Lorentz force law, with v = (electrostatics), contains more: the concept of a field E, the mediator 
of the force, so that action is not instantaneous. But there is one nice feature of Coulomb's law: it tells you that our 
space (not space-time) must be 3- dimensional. How? Try to write down in d-dimensions Gauss' law of constancy of 
electric flux over concentric spheres enclosing a point charge q, and you will get a Coulomb-like law: the field goes as 
l/r'*~^. That it goes as inverse square tells you that d = 3. 
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with 7 = (l — w^) . Note that Aq > 1 and deth. = 1; the transformations that satisfy these two 
conditions are known as proper Lorentz transformations. 

A note of caution here. One should write either A^,^ or Aj,^ to indicate clearly which one is the 

first index and which one is the second. Otherwise, one may write, for an antisymmetric mixed 
tensor, = — which is confusing to say the least. However, we will hardly fall in such a trap 
during this course, so can afford to be a little cavalier in the positioning of the indices. 

Eq. (7) defines a contravariant four-vector; the transformation law is 

A'^ = Ai^A" (9) 

which can also be written as {dx'^ ldx")A^ . Note that repeated indices are summed over. 

The invariant interval between two nearby space-time points in the Minkowski space is given 

by 

ds^ = dt^ - dx^ - dy'^ - dz^ = q^^dx^'dx'' (10) 

where r]^^^ = diag{\, —1, —1,-1) is called the metric tensor or simply the metric. People who deal 
with gravity and curved space-time call this flat space-time metric tensor and denote the full metric 
tensor by (7^^, but let us not go into that. 

Let's see what we get. 



• There are quantities with more than one Greek index (also called Lorentz index). Quantities 
with two such indices are called rank-2 tensors, r}^i, is the simplest example of that. One can 
construct higher rank tensors too. 

• The metric tensor can be used to lower the Lorentz index. We can define another four-vector 

= rjn^A'^ . This also has four components, but they are (B^, — B). So, the time component 
is unchanged (that's why whether you call it B'^ or Bq is absolutely immaterial) but the 
spatial components reverse sign ^. This is called a covariant four- vector (also known as a 
1-form) and transform as 

K = = ^5.- (11) 

Note that A'^ is just the inverse of A(^. 

• Eq. (10) also tells you something about forming a Lorentz scalar from two or more Lorentz 
vectors or tensors (vectors are nothing but tensors of rank 1). If in a product a Lorentz index 
occurs in pair, once as a superscript (contravariant) and once as a subscript (covariant) that 
particular index is to be summed over from to 3. This index becomes a dummy index; the 
process is called contraction. Thus, A.B = A^^B^ = AqBq — A.B is a Lorentz scalar. Note 
that A^^B^ = Afj_B^^. This is nothing but a dot product of two four- vectors. 

• Eq. (10) can be written as dx"^dx'^ = A.i^A.'^dx^dxa. If ds"^ is invariant, we must have A{^A^ = 
(5^. (Is it correct to write A(^Aj; = 1?) 

One can also form a contravariant metric tensor r]^^ to raise the index: 

= j^i^^A.. (12) 



®This statement depends on the choice of the metric tensor. Griffiths, for example, uses r)^v = diag{—l, 1, 1, 1), and 
in his case, only the zero-th component reverses sign. This is generally the choice that people who work with gravity 
make. The reason is that the curvature of a sphere in this metric comes out to be positive. On the other hand, the 
energy-momentum relation reads = — m^. Since I am a particle physicist, I will use the metric diag(l. —1, —1, —1) 
that keeps = m^, and will not bother about the sign of the curvature. For my case, timelike separations have 
ds^ > and spacelike separations have ds^ < 0, which is just the opposite of Griffiths. Anyway, before looking at a 
textbook, first check the metric it uses. 
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Since we can write ds^ = 'q^^^dx^dx^, and as x^ = (t, — x), it follows that 7]^^" is also 
diag{l, -1,-1,-1). Obviously, 'q'^'^'q^.a = 5^ and r]^''r]^„ = 4. 

One can form tensors of rank 2 (or higher) by taking products (not contractions) of Lorentz 
vectors; e.g., I can write A^^" = C^D". The transformation law is obvious: 

A'"" = A^A^A"^. (13) 

The transformation law of a covariant tensor of rank 2 is B^iy = A^A^i?^/?- There can be tensors 
of mixed type, with one (or more) contravariant and one (or more) covariant index; an example is 
A(^. You can easily formulate the transformation laws of higher rank tensors. 

A rank-2 tensor has 16 components, and is generally written as a 4 x 4 matrix. It can be 

symmetric (^^^ = A'^^^), antisymmetric {A^^'^ = —A'^^), or with no such obvious property. Metric 
tensor is symmetric; we will encounter the most important antisymmetric tensor in physics shortly. 

Examples of commonly used four-vector are (i) = (t,x), (ii) p'' = {E,p), (iii) = (p,j) 
(charge and current densities), (iv) A^^ = (i^. A) (scalar and vector potentials), etc. Use of four- 
vectors can make life more elegant; for example, the Einstein mass-energy relationship is just 

p^Pfj^ = rr?. (14) 

Another important four-vector is 9^ = d/dx^ = {d/dt,'V) [and = {d/dt, —V)]. Note that here 
the negative sign comes for the spatial part of the contravariant vector. (That's only natural: 
you expect df^x^ = 4, right?) The contraction of with any four-vector A>^ is known as the 
four- divergence ^. 

Q. Show that a symmetric rank-2 tensor has 10 independent components, and an antisymmetric 

rank-2 tensor has 6. 

Q. If A"'^ is a symmetric tensor and Ba/s is an antisymmetric one, show that A^^^B^p = 0. 

Q. How many components does the tensor A!^ have? How does it transform? 

Q. Show that the Lorentz transformation equations for the coordinates when the velocity v of the 

moving frame is in an arbitrary direction is given by 

T — 1 

t' = ^{t - v.x), x' = x -I- ^^(v.x)v - 7vi. (15) 

Q. Prom eq. (15), find the general form of A(^. 

Q. Show that A.B is a Lorentz scalar: A'.B' = A.B. 

Q. Suppose that, in a coordinate system, the metric tensor g^j^y is diag{l,—a^/{l — 
kr^),— a'^r'^,—a^r^ sin'^ 9) where a is a function of time. If g^vg^°' = what should be the 
form of 5^^? (This metric is called the Friedmann-Robertson- Walker metric and is used to describe 
the smooth expanding universe, but some texts use an overall minus sign in the metric.) 
Q. Consider the equation Q^f^ = SirGT"'^, where G is a constant and T"^ is a symmetric tensor 
satisfying daT°'l^ = 0. How many independent equations does this tensor equation represent? 
Q. All components of a four- vector must have the same dimension. How should you insert factors 
of c in p^ and in the SI system? 

Q. How should A^ and A, the components of the four-potential scale when you go from RLH to 
SI system? Show that A'^ = (A'^,cA) is the correct four- vector in the SI system. 

2.1 The Relativistic Action and Distribution Functions 

This is a digression and has almost nothing to do with the material that follows, but since we are 
going to discuss relativistic electrodynamics, the discussion will be incomplete without this. 

contravariant vector in the denominator of a fraction is equivalent to a covariant vector in the numerator for 
contraction. 
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In Newtonian mechanics, the equations of motion are obtained by minimising the action S = 
J Ldt. In rclativistic theories the volume element d'^x is invariant. This can be seen from the fact 
that the Lorentz transformation is a 'rotation' in the 4-dimensional space. Another way to see this 
is to check that the Jacobian is unity (do this). Anyway, I can write S = J Cd'^x where L = J Cd^x, 
£. being the Lagrangian density. If S is invariant, so is £. This formulation is more relevant for field 
theories. We will see later how to get the equations of motion (i.e.. Maxwell's equations) from the 
Lagrangian density of the electromagnetic field. In quantised versions, all particles are assumed to 
be excitations of their corresponding fields, so one only needs to know C there. 

For classical and rclativistic single-particle dynamics, the action must be only a function of s, 
where ds'^ = dt^ — dx"^ — dy^ — dz^ = dt^{l — v^) and hence ds = dty/l — v^- Let 

S = -a r ds = -a \/l - v'^dt, (16) 

Jb Jti 

where a is some constant. In the nonrelativistic limit (v <C 1) we have 

L = -aVl - «i -a + ^av^ H , (17) 

so by comparing with the nonrelativistic free Lagrangian (^mv"^) we see that a = m and hence 
S = —m J ds 

Now, ds^ = dsds = dx^dx^, so ds5{ds) = dx^5{dx^). To determine the dynamics, we vary the 
action and get 

/■\/ , ^ f'' dxi'Sidxa) „ ,ft /"^ dui" , 

Ob = —m / o{ds) = —m / — — = —m / u^doXfj_ = —mu'^ox^\g^ + m ox^——ds, (18) 

Ja Ja Ja Ja 

where u'* = dx^^/ds is the proper 4- velocity. If 5x^ vanishes at the end points, we obtain 

which is a generalisation of force- free motion in rclativistic mechanics ^. 

If we wish to deal not with a single particle but with a large collection of them, we introduce a 
distribution function in nonrelativistic mechanics (remember the Maxwell-Boltzmann distribution). 
We do the same for rclativistic motion. Let us consider a set of N particles, each of mass m, 
described by a distribution function f{p'^) at any given location in space. The total number of 
particles can be written as 

N= f d'^pB{pQ)5{p^ -m^)f{p^). (20) 



The delta function ensures that all particles are real and the theta function tells us that the energies 
are positive. Since A'', d^p, 9 and 8{p'^ — m^) are all Lorentz invariant, so is /. Now we can write 
the delta function as 

Sip" - m") = 6{j>l -E^) = ^ [5{Po -E) + 6{po + E)] , (21) 



^Thcrc is an extra term, — m, in the Lagrangian. But the Lagrangian is never unique and the equation of motion 
does not change by adding a constant term to L. However, it shifts H, the Hamiltonian, by +m: we now have the 
rest energy also apart from the standard kinetic energy. 

*The trajectory is called a geodesic. In flat space-time, it is a straight line, analogous to the ordinary force-free 
motion. If the space-time is curved, the metric depends on the coordinates, and the geodesic equation looks more 
complicated. It can be shown (and this is the central point of the General Theory of Relativity) that a locally non- 
inertial frame is equivalent to the presence of a gravitational field. Thus, gravitational effects can be taken care of 
by studying the force-free trajectory in a non-inertial frame. The reason for such an equivalence, or the equivalence 
between gravitational and inertial mass, is that the gravitational field equally affects all particles. One cannot make 
electromagnetism purely a property of the space-time because the charged particle trajectories depend on q/m and 
are different for different charged particles. 
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where E"^ = + m^. Because of the theta function, only the first delta function will be effective, 
and the po integration will just replace pQhj E, so 

N = \j ^f{po = E,p). (22) 

As N and / are Lorentz invariant, so is d^p/E. 

Finally, note that E = 7m and ds/dt = 7"^, so Eds/dt is Lorentz invariant. Multiplying the 
numerator and denominator by d^x, wc get Ed'^x{ds/d'^x). The quantity in parenthesis is Lorentz 
invariant, so Ed^x is also invariant. Taking the product, we see that the phase space element 
d^pd^'K is relativistically invariant, though none of the individual terms are. 

Q. Prove eq. (22). 

Q. Show that the Jacobian of Lorentz transformation is unity. 



3 Relativistic Formulation of Maxwell's Equations 
3.1 Field Tensor and Its Dual 

The electric and magnetic fields in eq. (5) can be derived from scalar and vector potentials: 

dA 

B = VxA, E = -V(/)-— . (23) 

They automatically satisfy V.B = and V x E = —d^/dt (Faraday's law), but the other two 
equations give interesting results. Ampere's law (with Maxwell's correction) gives 



while Gauss's law gives, after the addition and subtraction of d'^(f)/dt^, 

If you remember that = ~V^, you can easily combine eqs. (24) and (25) into 

d^d^'A'' - d'^d^A^ = f, (26) 

or in an even more elegant way: 

d^F^'''=f (27) 

where 

pui^ = di'A'' - d'^Ai' (28) 

is the electromagnetic field tensor, a rank-2 tensor which is, by construction, antisymmetric. This 
is by far the most important antisymmetric tensor in physics. If you have done the exercises in the 
last section, you know that this should have six independent components. Let's find them out and 
explicitly construct F^^^. 

The six independent components of F^^" are F^^, F^^, F^^, F^^, and F^^. For the first 
term 

^01 ^ ^0^1 _ 51^0 ^ (29) 
Similarly, = —Ey and = —Ez- From B = V x A, we have 



9 



Thus, 
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The covariant tensor Fj^i, is obtained by 
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E, 
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y 



(32) 



ie., by the substitution E — — E. The reason for this sign change is that we need one r/oo and one 
r}ii for the electric field components, whose product is —1. For the magnetic field components there 
are two —Is from the metric tensors, so there is no sign change. 

One can construct another rank-2 antisymmetric tensor from F^^^ . Define a rank-4 completely 
antisymmetric tensor e^naji as 



^ixvap = +1 for 0123 or even permutations, 
= —1 for odd permutations of 0123, 
= otherwise 



(33) 



and e^^°'^ = —e^ivap (this convention is opposite to that of Jackson, but the final results will be 
identical). Then 
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is called the dual tensor to F^'^ and is obtained by the substitution E ^ B, B ^ — E in eq. (31). 
Do you remember that we have talked about such a symmetry of the Maxwell's equations earlier? 
So what does C^" have to do with electrodynamics? 

It is quite straightforward to show that 



leads to the second pair of Maxwell's equations. 







(35) 



Q. Given eq. (31), show that eq. (27) gives two of the Maxwell's equations. 
Q. Show that eq. (35) indeed gives the other two equations of Maxwell. 

Q. As and G^"^ are both antisymmetric rank-2 tensors, one can construct three Lorentz scalars 
out of them: F'^'Fn^, G^'^Gn^, and F^'''G^y. Show that the first two are proportional to E^ - B^ 
while the third one is proportional to E.B. Hence these two quantities are invariant under Lorentz 
transformations. 

Q. Show that under the parity operation x — x, F^^^F^^ is invariant but F^'^G^y is not. 

Q. Why one does not construct a scalar like p^PuE'^'^ where p is some four-momentum (say that 

of the electromagnetic field)? 

Q. Convince yourself that €3012 = —1- What are the values of 62103, eo32i and £2032? 
Q. Show that rj^'^e^ua/s = 0. 

Q. Show that the transformation + d^X, where A is some scalar function of x, keeps F^'^ 

invariant. Such transformations are generically called gauge transformations; this particular one is 
called the Lorentz gauge. 
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3.2 Transformation of the Fields 



How does the electric and magnetic fields transform under Lorcntz transformations? It has been 
beautifully demonstrated in Griffiths how, from a physical perspective, one can understand the 
field transormations; in particular, if E (or B) is zero in a frame, how it can have a nonzero 
value in another frame. Thus, what appears as an electric phenomena in one frame may appear 
to be a magnetic phenomena in another frame (this led Einstein to his celebrated 1905 paper On 
the Electrodynamics of Moving Bodies that introduced STR). Thus, unless one always deals with 
nonrelativistic motions, it is better not to talk about E or B separately, but to talk about F^*^. 

We will derive the field transformation laws simply, from the transformation law of a rank-2 
tensor, eq. (13): 

= A° A^F°^. (36) 

Consider the motion to be along the common x-axis; thus V = (u,0,0). Only Ag, A?, Aj and Aj 
are nonzero, but = F^^ = 0, and F^^ = -F^^. This gives 

= A°A}F°i + A?AjFi° = (72 - 7^2) i^oi ^ ^37) 
Let us also work out a couple more: 

^0=^' = AOA^F"'^ 

= A'^AlF^^ + A'iAlF^^ = j{-E^-vBy). (38) 



F^^' = AlA^F"" 
^-B', = AlAlF"^ + AlAlF'^ = jivEy-B,). (39) 

The complete set looks like 

E', = E,, E'y = j{Ey-vB,), Ei = j{E, + vBy), 

B'^ = Bx, B'y = j{By + vE,), B', = ^{B,-vEy). (40) 

Eq. (40) is symmetric under the exchange of E — B, B — — E. So you will get the same 
transformation laws if you start from the dual tensor G^^. 

Q. Show that for a general Lorentz transformation 

2 2 
E' = 7 (E + V X B) - ;^;^v (v.E) , B' = 7 (B - v x E) - (v-B) . (41) 



3.3 Fields due to an Uniformly Moving Particle 

Eq. (40) tells us how to get the fields in a comoving frame. Let us now look at the other side of the 
coin: suppose a particle is moving uniformly in the lab frame (say along the positive x direction); 
what should the fields be due to that particle? 

The approach is very easy; go to the frame where the particle is at rest, calculate the fields 
there, and boost back to the lab frame. What we need is the reverse transformations of eq. (40). 

Suppose we are in the frame 5', and the particle with charge q moves with a velocity v along the 
X-axis. The detector is located at the point (0, b, 0) in the S frame (this can always be arranged; 
call the nearest distance of approach b and let the y-axis pass through the detector). Let t = t' = 
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Figure 1: Fields in S due to a uniformly moving particle which is at rest in S'. 



when the particle was at the origin of S, i.e., at the nearest distance of approach. Let S' be the 
frame where the particle is at rest, so the two frames coincide at t = t' = 0. Also, let n be the unit 
vector along the line joining the instantaneous position of the charge (at the origin of S') and the 
detector, and n.v = cos^^. Thus, 6 = rsin^ and vt = —rcosip. 

At a time t in 5 and t' in S', the coordinate of the detector in S' is x[ = —vt', x'2 = b, and 
X3 = 0. The distance is r' = ^ [vt'^ + The electric field components in S' are 

qvt' I qb 



E[ = - 



Ei = 0. 



(42) 



The magnetic field components B[,B'2,B'^ are all zero. 

In unprimed coordinates, the fields in S' (I have not yet implemented the field transformations) 
look like 

E[ = - . ... ^^f K = . ... * = 0. (43) 



47r (62 +^2^2^2)3/2' 2 + ^2^2^2)3/2 ' 3 

Now we boost the fields back to the lab frame with the inverse of eq. (40), i.e., replacing v by —v: 

q^vt 



B: 



Ei = E[ 
E-2 = 7-E'2 
jvE'2 



47r (62+^2^2^2)3/2' 

jqb 

47r (62 + 72^;2t2)3/2 ' 

^vqb 

47r (62 + 72^;2t2) 



3/2 ■ 



(44) 



This we can write in a more compact way. Note that E1/E2 = —vt/b, so E is always directed 
along n, the unit radial vector joining the present position of the charge to the detector, just as a 
static Coulomb field. Also, the denominator (62+72^72^2^3/2 written as r^7^(l — ^2 sin2 



(replace 6 and vt by r and ip and use 7 — 1 



E 



72t|2)^ SO 

qr 



and 



47rr^72 (1 — ^2 sin2 ij;^ 



B = V X E. 



3/2 



(45) 



(46) 



Thus the passing charge induces a magnetic field. That's nothing new; in fact, for low velocities 
(7 1) the expressions for E and B are identical to the ones that you get from the Coulomb 
and the Biot-Savart law. More interesting is the high-velocity limit {v 1,7 » 1). At t = 0, 
transverse electric and magnetic fields E2 and S3 show a sharp peak, of approximate height jq/b'^, 
and the peak falls off sharply as even for a small t the 62 term in the denominator will be utterly 
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negligible. This is indistinguishable from the effect of a pulse of plane polarised radiation moving 
in the x-direction. The longitudinal field Ei varies rapidly from negative to positive, has a zero 
time integral, and is actually zero at t = 0; the detector won't feel its presence. 

Let us note an important point right here. If you calculate the Poynting vector E x B from 
eqs. (45) and (46), you will get a nonzero result, so the field carries some energy. However, it does 
not radiate. We will see later that radiation means the presence of energy at an infinite distance 
(as if the energy decouples itself from the charge configuration and moves to infinity). Here if you 
integrate the Poynting vector over a surface at infinity, the answer will be zero, since both the fields 
fall off as (and hence the Poynting vector as 1/r^) while the surface area of the sphere only 
grows as r^. Thus, to get a radiation, the fields cannot fall off faster than 1/r, and you will soon 
see that only if the charge is accelerated, there is a 1/r component in the field. Thus, a charge 
moving with constant velocity never radiates; to get radiation one must have an accelerated charge. 

There is another very elegant way to see that a charge moving with uniform velocity won't 
radiate. A static charge does not radiate; E oc and B = 0. But a charge moving with uniform 
velocity can be made static in another inertial frame moving with the same velocity. Physical laws 
must be invariant in all inertial frames, so such charges do not radiate. There is no cheating here; 
this is as watertight an argument (if not better) as the one given earlier, or the one that we will 
see later using retarded potentials. // you know the full power of the Special Theory of Relativity, 
a lot of relations can be deduced very quickly. 

Q. A charge e moves along the x-axis with a constant velocity v and passes the origin at t = 0. 
How do you write the form of j? [Hint: you need some Dirac S-functions.] Show that for v <C 1 
you indeed get the Biot-Savart law. 

Q. Show that the full width at half maximum of the transverse pulse is of the order of b/jv. 
Q. What should be the form of eq. (44) in the SI system? 

4 Lorentz Force Equation and Its Generalisation 

We know that x is the spatial component of a four-vector x^^. Is dx/dt the spatial component of a 
four-vector? The answer is no, since both x and t change under Lorentz transformation. Rather, 
we define a proper time r which is the time that passes in a clock of the system where the particle 
is at rest. Clearly, this is a unique quantity and does not Tindcrgo any Lorentz transformation. 
So the hybrid quantity dx/dr behaves like the spatial component of a four- vector ^. This is the 
four-velocity, which we will call u^, and whose zero-th component is dt/dr = 7. On the other hand, 
V = d'x./dt is the ordinary three- velocity. Remember that it is wrong to talk about an ordinary 
four- velocity v^. 

Similarly, ordinary force F = dp/ dt is not the component of a four- vector. Rather, K = dp/dr 
is. Obviously, K = 7F and = dE/dr is the power (call it proper power if you like). 

Lorentz force is an ordinary force. Let us find the four-force K^^ (sometimes called the Minkowski 
force) which behaves as a Lorentz vector. The force should be constructed from the products of 
fields and velocity components, and the only option is Ui^F'^". There should be the electric charge, 
so we throw q also, and try with the product qu^F^'^ . 

Wait. First, why don't we start with the dual tensor G^^'^l Well, we wish to get the Lorentz force 
in the nonrelativistic limit, but if we started from G^'^, we would have ended up with g(B — v x E), 
something which we don't want. Second, is the electric charge Lorentz invariant? The answer is yes; 
experimentally we never saw a fast-moving electron to have any charge but —1.6 x 10~^^ coulomb, 
and theoretically there arc good reasons why it should be a Lorentz scalar (there are effects that 
change the electric charge, but that has nothing to do with Lorentz invariance). More about this 

®I call it hybrid since x is measured in S and r in S". 
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later. 

The proper velocity is given by u^, = (7, — 7v), and the spatial components of can be easily 
computed: 

Ki = q (uqF^^ + U2F^^ + usF^^) = 97 (E + v X B)^ . (47) 

Thus the Lorentz force law is restored at the nonrelativistic limit 7 1. Kq = q'yv.E gives the 
rate of change of energy of the particle. 

Another way to get the equation of motion in an electromagetic field, analogous to that of eq. 
(19), is this. The action can always be written as S = —m J C{x^, u^)ds. We can expand C (which, 
to maintain translational invariance, should not better be an explicit function of x^) as 

C{x^', n^) = 1 + CQct){x) + ciA^{x)u^ + C2g^^{x)u^v!' + • • • , (48) 

where the argument x of the fields is a shorthand for the four- vector x'^ . 

In practical applications, terms upto C2 are enough. The curved space-time metric tensor g^^, 
is important for gravitational interaction. The scalar term can also be included in this by adding 
Vfj.i^'t'ix) to gfj,u, as u^^Ufj, = 1. We are left with the ci term, which is relevant for electrodynamics. 

Let us write the action as ^ 

S= f {~mds - qA^dxf") , (49) 

J a 

where we have substituted q, the charge of the particle, for ci since ultimately we must have the 
Lorentz force law. This tells us the total Lagrangian is L = — mVl — v'^ — qA'^ +qA.v as dx./dt = v. 
Varying the action with respect to the trajectory x'^{s), we get 

SS = - (""^^^^^ + + qSA^'dXf,^ = 0. (50) 

Using SA^ = {dA^/dx'')6x'' = d^A^Sx'', and dA^ = d^A^^dx", and A^diSx^") = d{A^6xi^) - 
dj^A^dx'^Sx'^, we get 

ds - {mu^ + qA,,) 6x^\^^ = 0. (51) 

In the third term, we interchange the dummy indices n and z/, which changes nothing. The variations 
for which Sx'^ vanishes at the end points lead to the equation of motion: 

= qF^,u\ (52) 







m- 



dUf, 
ds 



fe^ + q{d^A^)u''6x'' - q{d^A^)u^'6x' 



4.1 Motion in Combined Uniform and Static Electric and Magnetic Fields 

Consider a particle moving in constant uniform electric and magnetic fields. They need not be 
parallel; in fact, here we will only discuss the case where they are perpendicular. In such a combi- 
nation, K'^ is not zero, so the particle receives energy from the field and its velocity always changes. 
Still, we will try to make life simpler by going to a frame where one of E and B is zero. Since 
E.B is a Lorentz invariant quantity, if the fields are not perpendicular, the dot product will always 
retain the same value in different frames and it is impossible to make one of them zero. Moreover, 
since E^ — B^ is also Lorentz invariant, only that field which is smaller in the starting frame can 
be made zero. 

Suppose |E| < |B|. Also suppose E = (0, 0) and B = (0,0,5). Let us find the components 
of the fields in the frame S' where E is zero. The velocity of this frame should be in the x-direction 
(this keeps and zero); actually, one needs to choose v as 

(53) 
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Here v = {E/B, 0,0) and 7 = B/^B^ — E'^. Note that |v| < 1, consistent with the second 
postulate of STR. With such a choice, the fields in 5" are 



Thus, E vanishes in S' and B gets scaled down by a factor 7 while retaining the same direction. 
The motion in S' is trivial: a circular motion around the lines of force. The direction in which the 
spiral turns depends on the charge of the particle. When we come back to the frame 5, the linear 
motion along the positive x-axis is superimposed on the circular motion (this is called an E x B 
drift), so we get a true spiralling trajectory. 

The particle always moves in the positive x direction, regardless of its charge, since its charge 
entered nowhere in the calculation. The only difference will be the direction of its turn, whether it 
is a right-handed screw or a left-handed one. 

Such crossed fields can be effectively used to get monoenergetic particle beams. Suppose we get 
a number of electrons from a cathode-ray tube, with a velocity distribution. If we pass them through 
such crossed static fields, only those electrons with velocity equal to E/B will go undeflected; others 
will bend towards or away from the electric field. If we know what energy is needed, we can set 
E and B accordingly. Suitable entrance and exit slits with momentum selectors like deflecting 
magnets can produce a very pure monoenergetic beam. This is particularly important in the high 
and low-energy collider experiments. 

If |E| > |B|, there is no such spiralling motion. It is easy to show that with a velocity v = 
(E X B)/£'^ one can make B = in S" while E' will be E/7. This is the motion of a particle in 
a constant electric field; it gets continuously accelerated. In S the overall velocity of S' will be 
superimposed on it. 

Q. Get the components as given in eq. (54). Show that E.B and E^ — B^ are indeed Lorentz 
invariants. 

Q. If E and B are perpendicular in a frame, show that they will remain perpendicular in all co- 

moving frames. 

Q. If E and B make an acute angle with each other in a frame, show that the angle between them 
will always remain acute. 

Q. What should be the form of v (eq. (53)) in the SI system? 



5 Lagrangian and Equation of Motion 

We can write the Maxwell's equations in terms of E and B, or in terms of F^^'^ and its dual, or 
in terms of the four-potential A^^ = (^, A). All formulations are completely equivalent. The first 
one uses the quantities you are familiar with, the second one is elegant and manifestly Lorentz 
covariant, and the third one (which is just a longhanded way of the second; Fi^" = d^A^ — d'^ A^) is 
suitable to get the equations of motion. What should the equations of motion be? Well, you know 
the answer: the Maxwell's equations! But we will try to get them from a different starting point. 

Prom our lessons in classical mechanics we know that the Lagrangian L = T — V \s a function 
of the generalised position q and the generalised velocity q. It can even be an explicit function of 
the time t, but let us consider only systems where L = L{q,q). The action S, defined as J L dt 
between the initial time ti and the final time t2, must be an extremum for an allowed path. Prom 
this follows the variational principle 



E'^ = E^ = 0, E'y = j{E - uB) = 0, E'^ = 0, 
B'^ = B^ = 0, B'y = 0, B', = ^{B-uE) = -. 



(54) 




Ldt = 



(55) 
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and the equation of motion through the Euler-Lagrange prescription: 



1(1)- 1 

The Hamiltonian is H{q,p) = pq — L where p = dL/dq. 

What happens for a classical field which can theoretically be extended over an infinite volume? 
There is a chance that the quantities that we wish to evaluate may turn out to be infinite; in fact, 
the total energy of the electromagnetic field, / ^(E^ + B^)dt; is infinite if we consider an infinite 
volume with a constant electric or magnetic field. But this is not a defect of the fields; the result is 
infinite just because we integrated over an infinite volume! Thus it is better to talk about a density, 
e.g., the energy density 6, where the total energy is / £^ dv. Evidently, S is finite. 

We can similarly talk about a Lagrangian density C, with 



/ 



Cdv = L. (57) 



Apart from C being finite, this has an extra advantage; the action looks better from a relativistic 
point of view: 

S = j d'^x. (58) 

What should £ be a function of? The generalised coordinate is the field tp (do not confuse it with 
the scalar potential 4>) which in turn depends on the coordinates x^, and the spatial and temporal 
derivatives of ip, namely, d^tp (again, from a relativistic standpoint, if we introduce a temporal 
derivative we should at the same time introduce the spatial derivatives too). The Euler-Lagrange 
equation is slightly more complicated, but analogous to eq. (56): 

d fdC\ ^ f dC \ dC ^ 

- I — 1 + V. ( ^,7^ ] - _ = 0. (59) 



dt\dp) \d{Vp)) d^p 
This can be expressed in a more compact notation: 

o dC dC ^ 

'^d{di^p) dp 

This gives the equation of motion of a field. 

Before we go to the dynamics of the electromagnetic field, let us get some practice with some- 
thing far easier: the motion of a charged particle in an electromagnetic field. 



5.1 Charged Particle in an Electromagnetic Field: the Generalised Momentum 

The ordinary 3-dimensional Lagrangian must be a scalar; we don't want it to change under some 
coordinate transformation. Similarly, a relativistic Lagrangian (we will talk about the Lagrangian 
density only for the fields, not for individual particles; this is nothing but classical mechanics with 
the force provided by the electromagnetic fields) must be a Lorentz scalar. I don't know what the 
Lagrangian is (Lagrangians are, as a rule, never derivable from first principles, they are not even 
unique; their only test is the reproduction of the equation of motion) but in the nonrelativistic 
limit it should look like —ecj) (I am using e as the charge of the particle, to avoid confusion with 
the generalised coordinate q). But 4> is the zero-th component of A^^, so the actual 4-dimensional 
Lagrangian must involve A^^. What should it be contracted with? I have only two four- vectors, x^ 
and the proper velocity u^; but the Lagrangian must be translationally invariant too, apart from 
being Lorentz invariant, so it cannot contain x'^. Thus, the interaction Lagrangian must be 

Lint = - -u^A^ = -e0 + ev. A, (61) 
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where the factor of 7 has been introduced to get the nonrelativistic Lagrangian as — and not as 
—67^; this factor also changes u to the ordinary velocity v. 

Eq. (61) modifies the canonical 3-momentum P from the mechanical 3-momentum p by the 
addition of a term eA (remember that P = dL/dv): 

P = p + eA. (62) 



The total energy is the mechanical energy, vp^ + rr?, plus the electrostatic energy ecp: 

W = y(P^^eA)2Tm2 + 6(1). (63) 

This is nothing but the usual energy-momentum relation p'^p^ = vr? with 

p = [E, p) = (M^ - e(/), P - eA) . (64) 

Without the electromagnetic field, p = P and = so the effect of the field is to make the 
replacement 

p'' ^p^'- eA^. (65) 

Eq. (65) is the most important equation in the theory of interaction of electromagnetic field with a 
charged particle. Making this minimal substitution gives the entire dynamics; of course, one has 
to rewrite A^^ in terms of E and B. 

In fact, eq. (65) has a much greater impact on physics. You know that quantum mechanics 

stems from classical mechanics by the identification of the operator — iV with the 3-momentum p 
and the identification of the operator id/dt with the Hamiltonian or the energy E (we use ?i = 1); 
in short, by the identification of id^ with p^. For an electron in an electromagnetic field, substitute 
id^^ with id^ + eA'^ (electron's charge is negative) and operate this operator on the electron wave 
function, and you get the entire quantum electrodynamicsl 

Q. We have not talked about the gauge invariance of the Lagrangian, but it is imperative to check 
that, particularly when the Lagrangian is not a function of E and B but of which is not a gauge 
invariant quantity. First, check that if we add a total derivative df{x,t)/dt to the Lagrangian the 
equation of motion remains unchanged. 

Q. Consider a gauge transformation Af^ ^ A^^ + d'^X{x). Show that eq. (61) is changed by a total 
time derivative. You will require 

so better be convinced that you know this! Anyway, the Lagrangian is not invariant, but the equa- 
tion of motion is. 



5.2 Lagrangian for the Electromagnetic Field 

As we have mentioned earlier, in this case we will consider the Lagrangian density C and not L. 
The Lagrangian density must be a Lorentz scalar, because we do not want it to get transformed 
under a Lorentz transformation. This follows exactly the same logic that in ordinary 3-dimensional 
case, the Lagrangian must be a scalar and can never be a vector. Moreover, electromagnetism 
respects parity. This means that if you notice the motion of a particle in an electromagnetic field, 
and then reverse all the spatial coordinates of the system, the new motion will also be allowed. 
Thus, we expect C to be invariant under parity transformation x — x too. 

These two considerations severely limit the possible options. If we consider the Lagrangian 
density for a free field, j^^ = 0, and we have no other quantities at our disposal except Ff^" and 
G^^'^. Thus, there are only three possible choices: F^^'^Ff^^,, C^^G^i,, and F^^'^Gf^^- Among them, the 
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third term, which is proportional to E.B, is not invariant under parity, since B is an axial vector 
while E is a proper vector. The first two terms are equal, but they satisfy both the criteria for 
being a valid term in the Lagrangian density. Thus, L oc F^^F^j^u. 

Let us start with 

C = —F^^F^,. (67) 

The factor of —1/4 is a matter of convention; any multiplicative factor would have given the same 
equations of motion. Using the explicit form of F'^'^ , this becomes 

£ = _ 1 {di'A^d^A, - d^^A'd^A^) . (68) 

We treat A^ as the electromagnetic field, not E or B. This has some justification: ultimately we will 
quantise A^^ and get the photon as the excitation quantum of the field. Anyway, eq. (68) involves 
only derivatives of A^^, not A*^ itself, so eq. (60) becomes 

Let us now compute dC/d{d^Au). The first term on the right-hand side of eq. (68) is just the 
square of d^A'^, so this will yield —d^A"^ (be careful about the position of the indices). For the 
second term, we note that 



d{dpAr) 



= \5i8lrr\^-d^A^ + \d^A,rr\^-bi6l 



so that 



\^^P-rf^dy^A, + l^d^A^yTPrf^ = A" (70) 



= -di^A" + d^Ai" = -F"" (71) 



d{d^A,) 

and the free-field Euler-Lagrange equations become 

di.Ft"' = 0. (72) 

This set of four equations are nothing but two of the Maxwell's equations, Gauss' law and 3- 
component Ampere's law, written in the absence of any external charge or current densities. Ob- 
viously, if we started with 

£ = -1g^-G^, (73) 

we would have obtained the other two equations of Maxwell. 

Why, then, we write the Lagrangian density in terms of and not its dual? Suppose we have 
a nonzero external four-current density = (p, j). In this case, we can write another term C, of 
the form j^A^, and 

jr = AF^-F^^-fA^ (74) 
yields the correct equations of motion, namely 

d^F''''=r. (75) 

There is no such magnetic analogue of j'^. The factor of —1/4 is needed to get eq. (75). 

Before we end this subsection, let us comment on the gauge invariance of eq. (75) (gauge 
invariance of electromagnetism is such a sacred principle that we want to check it at all steps). 
Consider the transformation Ai^ ^ A'^ + S'^A. Ff^" is invariant by construction and so there is no 
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problem with eq. (67). In eq. (74) the term j'^An is apparently not invariant, but gets an extra 
contribution of j^d^X. However, 

fd^X = d,{fX) - {d^f)X. (76) 

The first term is a four-divergence and vanishes when we compute S = J Cd'^x. By an analogue 
of the 3-dimensional divergence theorem, we can reduce the four-divergence integral to / j^X dv; 
but nothing goes out of the total volume, so the contribution is zero. The second term is zero 
only because the electric four-current is conserved by the continuity equation: S^j'^ = 0. Thus, 
the Lagrangian density of the electromagnetic field in the presence of an external current is gauge- 
invariant only because the current is conserved! The argument can of course be turned the other 
way around: the current is conserved because we demand gauge invariance. In short. 

Current conservation •s-^ Gauge invariance. (77) 

This is such an important statement that we give it a separate equation number. 

Why don't wc never talk about a term like A^^A^l This would have been a perfectly valid term 
in vC, but unfortunately this does not respect gauge invariance. For the electromagnetic field this 
does not matter, since such a term would give rise to the mass of the field quantum, and we know 
that as far as experimental accuracy goes, the photon is indeed massless, so there is no harm if 
the theory cannot accomodate a mass term. But there are theories where we need to have massive 
photon-like objects without breaking the gauge invariance, and a consistent formulation of such 
theories is indeed subtle. Unfortunately we don't have time and space to discuss that here. 

Q. Show that A^A^ is not gauge invariant. 

Q. Prom eq. (74), show that the momentum conjugate to A^ is zero. 

Q. Get the momentum conjugate to A. For this, compute dC/d{cfA^^). Show that this is equal to 

the electric field E. 

Q. If we start from C = —\G^^G^i,, will you expect the momentum conjugate to A be the magnetic 
field B? Explain. 

Q. If A^ gives rise to the photon field, we expect the degrees of freedom of A^ and a real photon 
to be equal. Are they equal? 



5.3 Energy and Momentum of the Electromagnetic Field: Poynting's Theorem 

Consider a free electromagnetic field {i.e., = 0). Let us construct a rank-2 contravariant tensor 
T"^ as 

where C is given by eq. (67). From our knowledge of classical mechanics, T'^^ is something like the 
Hamiltonian density of the field. Using eq. (71), we recast eq. (78) as 

rpap ^ -F'^^'^Pd^Ax - rf^H. (79) 

T"'^ is called the canonical stress tensor. If you remember (if not, check) that C = |(E^ — B^), it 
is easy to calculate T^^: 

-F^%Ax B^) = -E.aoA - ^ (e^ - b2) . (80) 

However, E = - - OqA, so -E.^qA = E^ -f E.V^. Putting this, 

= ^ (e^ + B^) + V.(<^E), (81) 
19 
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where we have used the free-field equation V.E = to get the last term on the right-hand side. 
Similarly, 

= -E^d'Aj = E^diAj. (82) 

Again, we note that 

(E X B) . = €kij€kimEjdiAm = {SuSjm - SimSji) EjdiAm = {EjdiAj - EjdjAi) , (83) 

where we have been a little cavalier in the positioning of the indices, but that does not matter since 
they are all ordinary three- vectors. Now we can write eq. (82) as 



-lOi 



(E X B) . + V. (AiE) . (84) 



Integrating eqs. (81) and (84) over the whole volume, the three-divergence terms drop out due to 
the divergence theorem, and we have 

J d^x = ^J (e2 + b2) d^x = E, J T^' d^x = J (E x B) d^x = Pi. (85) 

These are the usual expressions of the total energy and momentum of the electromagnetic field. 

From cq. (85), it is easy to formulate a differential conservation law, a four-dimensional analogue 
of Poynting's theorem: 

OaT"/^ = 0. (86) 
This is easy to prove. From the definition of T"^, 



dA^ - d{dM 

= df^£{Ax,d'^Ax)-df^£ = 0. (87) 

In the intermediate steps, we have used the Euler-Lagrange equation, the chain rule of differentia- 
tion, and the fact that JC is only a function of Ax and its derivative. 

However, note that T°'^ is not symmetric. It is also neither traceless (T"" ^ 0) nor gauge 
invariant. The traceless property is needed for a massless photon to emerge after quantisation; the 
symmetric nature is needed to conserve angular momentum of the field, and of course we would 
prefer quantities to be gauge invariant. To this end, we write 

= (^F" V"^/.:^ + ^^/"''-^^'''i^/.z.) - i^"''/'^^^, • (88) 

The quantity in parenthesis is symmetric. The last term can be written, with source-free Maxwell's 
equation ^|^F"^' = 0, as (f^'^A'^) . This is a four-divergence and hence gives zero under integra- 
tion. So we may neglect this term and define a symmetric stress tensor 6"''' as 

Qa/3 ^ (^F^t^^'^^F^, + Iry^'^F^-F^,^ . (89) 

Q. Satisfy yourself that you understand all steps that led to eqs. (81), (84), and (87). 
Q. T°^^ is not a symmetric tensor. Calculate and show that it is different from T^. 
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Q. Show that da@'^^ = 0. 
Q. Show that 



1 



2(E' + b2), 



@oo 

QOi ^ Qio = (E X B) . , 



EiEj + BiBj — -(5ij(E^ + B^) 



(90) 



Identify Q^^ with the negative of Maxwell's stress tensor. 



6 Potential Formulation 

The four-potential is not unique; the source-free Lagrangian density is invariant under a gauge 
transformation, and even the source term is if the current is conserved (and it is conserved) . 

If we want a unique solution for A'*, we must specify the gauge, and in relativistic electrody- 
namics, the Lorentz gauge, 5^^^^ = 0, is the most convenient one to use. In this gauge, the equation 
of motion is 

□A^=j^ (91) 

This is called the inhomogeneous wave equation. In the absence of any source, = 0, the equation 
reduces to the homogeneous wave equation. 

The question we would like to ask: what is the potential at any point due to a moving charge? 
Now that the gauge is specified, this is indeed a meaningful question with unambiguous answer. 

Is our counting of the number of degrees of freedom correct? Apparently, A^, being a Lorentz 
vector, has four components, while a real photon has two (left and right circular polarisation). 
Obviously, two degrees of freedom are just spurious; you may call it an artifact of our formalism. 
The gauge condition gives one constraint and hence removes one degree of freedom, but what about 
the second one? Well, if you look at the Lagrangian density of the field (even when a source term 
is present) you will immediately see that the momentum conjugate to A^ is zero, since there is no 
doA'^ term there. So the scalar potential is a cyclic coordinate and can be removed by equations of 
motion; it does not represent a true degree of freedom. 



6.1 Retarded Potential 

This is really a very simple concept: electromagnetic wave, i.e., light, takes a certain amount of 
time to reach the observer from the source. Thus, what the observer sees right now (say, at t) is 
some configuration that was there at some earlier time (say, at to). Obviously, to < t; in fact, there 
is no way to know what is happening to the source of the electromagnetic wave (let us just call it 
the electron) at the present moment t. 

If you happen to be familiar with the concept of the light cone, you know that points which are 
inside the light cone can be causally connected (i.e., there is a cause, the signal transmits, and the 
effect follows). If the signal travels at the velocity of light, it is the boundary of the light cone that 
gives causally connected points. If I draw the world line of the electron, only at the point where it 
intersects the light cone of the observer can I have any information about it. To be more specific, 
if I start my light cone from (to, xo) and it intersects the world line of the electron at {t, x), then at 
the latter point I can have the information about the electron when it was at the position {to, xo); 
this is the retarded position (and the retarded time) of the electron. Since the separation between 
these two points is light-like, 

R''R^c = 0, (92) 
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where = {t — to,x — xq). The point of intersection is, of course, unique. It is only the retarded 
potential, A^{to,xo), that we can talk about and calculate. 

There is another point of intersection on the backward light cone, but that is not of any physical 

relevance, since that violates causality; effect precedes cause. These are called advanced time and 
advanced position of the electron. However, they are not just esoteric concepts; it may be shown 
that when we quantise the field, the advanced coordinate is related with the motion of antiparticles. 

Q. Draw the light-cone diagram that we have discussed just now. Convince yourself that causality 
is satisfied only on the forward light cone. 

Q. How should the light-cone condition, R^Rfi = 0, look like if the velocity of the signal is some 
V <1 (remember, we take c = 1)? 



6.2 Digression: Lorentz Invariance of Electric Charge 

Why the electric charge e is Lorentz invariant? It would have been a total mess, both theoretically 
and experimentally, if it were not; but that is hardly an answer. Note that the charge density p is 
not Lorentz invariant, since it is the zero-th component of the four- vector j^. In a frame where the 
configuration is static (we neglect the internal motion of the configuration), j^^ = {f)Q,{)), and in a 
frame which is moving with a uniform velocity v with respect to the former, j'* = (p, pv); current 
is nothing but motion of charge. But j^j^ is Lorentz invariant, so 

P = 7Po- (93) 

To get the total charge I have to integrate over the volume, but the volume element dVo gets 
contracted to dV/^ (contraction occurs only in the direction of motion), so / pdV is a constant. 



6.3 The Lienard-Wiechert Potential 

Consider an electron of charge e (here e is negative) to be at a point Xj*, and the observer at xl^. At 
the proper frame, where the electron is at rest, we know the solution of the inhomogeneous wave 
equation: 



where ro is the spatial distance between the charge and the observer. This is nothing but Coulomb's 
law, but the retardation condition tells us that 

i?'^ = (x2-a;i)'^ = (ro,ro). (95) 

Our task is to write the solution of A^^ in some Lorentz covariant form which in the static limit 
reduces to eq. (94). Since the proper velocity is (1,0) in the static frame, so that u.R = ro, we 
can write 

At: u.R' ^ ' 

subject to the condition R? = 0. In a moving frame, = (7, 7v), and R^ = {r, r), so 

A'' = -rT-^ r(l,v). (97) 

47r(r - r.v) 

Eq. (97) is the Lienard-Wiechert potential. 

This is the simplest way to derive the form of the Lienard-Wiechert potential, though a clever 
guess is required. There are two other ways: one is the so-called method of 'information-collecting 
sphere', which introduces the concept of the retarded time in a roundabout way and does not 



22 



contain an iota of more physics insight than this (see Panofsky and PhiUips, or Griffiths, for this 
method). The second method is through the use of Green functions (see Jackson for a detailed 
analysis). This is no doubt more rigorous, but at the same time mathematically more complicated. 

Q. What should be the form of eq. (97) in the SI and the Gaussian systems? 

6.4 Fields due to a Moving Charge 

Suppose a charge is moving along some trajectory. What should be the electric and the magnetic 
fields at a certain distance? Well, to the observer, the present position is not known, so the fields 
should come out as a function of the retarded position only. We have to use the usual definitions 
of E and B, namely, E = —VA^ — dA/dt and B = V x A, but we have to use eq. (97) for the 
potentials. The extra complication comes from the r.v term in the denominator. 

There is one exception: when the charge moves with uniform velocity. Only in this case, if we 
know the retarded position and the velocity of the charge, the present position is also known, and 
we can express the fields in terms of its present position. But this is precisely the same game that 
we have played earlier! We have seen, from eqs. (45) and (46), that the fields are indeed directed 
towards the present position of the charge Let us recall the main conclusions: 

• A moving charge induces a magnetic field; 

• For high velocity the field is like that of a plane transverse wave; 

• The field does not radiate. 

Q. Explicitly show how you get eq. (45) from eq. (97). If you cannot, see Panofsky-Phillips, section 
19.2. 

6.5 The Fate of the Potential 

I have said enough to confuse you thoroughly. First, I said that the potential smells of action- 
at-a-distance and hence can be a valid concept only in a nonrelativistic theory. This is true for 
both scalar and vector potentials. When we bring in the Special Theory of Relativity, we also say 
that nothing can move faster than light in vacuum, so no signal of electromagnetic disturbance can 
travel to a far observer instantaneously. This kills such action-at-a-distance theories. The electric 
and magnetic fields are the only measurable and relevant quantities. Indeed, we can write F^'^ in 
terms of E and B alone. We can calculate the radiated energy in terms of E and B. Where, then, 
is the place of the potential? Why should we waste so much time talking about Why should 
we derive the form of the Lienard-Wiechert potential? Isn't it a useless concept in relativistic 
electrodynamics? 

To answer this apparent paradox, let us first understand what potential is. It is just a function 
of space and time. Its main use is that it gives the correct E and B, when a particular prescription is 
followed. It is not unique; electromagnetic gauge transformation tells us that. In the nonrelativistic 
case, the introduction of the potential helps us to calculate the scattering amplitudes (we have not 
shown how to calculate them in classical relativistic electrodynamics, and this is not an easy job). 
That's the end of the story. What we derive in the preceeding subsections is a function of space 
and time from which F^^, and hence E and B, may be calculated. 

But if that were the only case, we would scarcely spend so much time in deriving just another 
function. The fact is that potential — by now I mean A^^, the four-potential — comes back with 

^°Be careful about the notation. In eqs. (45) and (46), r denotes the present position of the charge. In eq. (97), r 
denotes the retarded position. They are obviously not the same. 
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a vengeance in quantum electrodynamics. This is a quantum field theory which describes the 
interaction of electrons and photons. In quantum mechanics, -0 is a wavcfunction that describes a 
particle, say an electron. In quantum field theory, 'ijj is upgraded to the status of an operator. This 
operator, acting on a state, may create an electron (so that the number of electrons is increased by 
one). It can also destroy a positron, but that is a separate issue. In short, ip is an operator — this 
we will call the electron field, and this is a function of x*^ — whose excitations are the particles. 

A^^ is the field of the photon; you quantise it and get the photons as field excitations. That's 
where its importance lies. You can still call it a four-potential if you like, but remember what 
Shakespeare said about roses. 



7 Accelerated Charge 

In this section, to avoid confusion, we will use the following convention: 

r denotes the retarded position of the moving charge, and yq denotes its present position if it 
would have moved with uniform velocity. Since it docs not, the present position can be something 
completely different, but that really does not matter, since any information about the present 
position is inaccessible to us. The velocity and the acceleration of the charge is v and a respectively. 
In general, p denotes the magnitude of the three-vector p. The angle between r and v is 9, while 
that between tq and v is ■0. 




Present pos. 



Retarded pos. 



Observation pt. 



X 



Figure 2: Retarded and (virtual) present positions of a particle. 



We call ro the virtual present position of the charge, and clearly fq = r — rv (remember that 
V < 1 and light moves with unit velocity). The Lienard-Wiechert denominator is 



r.v 



ro'\J 1 — v'^ sin^ '0, 



(98) 



(To get this, use r^ = r"^ + r^v^ — 2rr.v and = + (r.v)^ — 2rr.v, and subtract one from the 
other; also use rxv = roXv.) The expressions for E and B for a uniformly moving charge is 
given in eqs. (45) and (46). In fact, we would like to write 



B = vxE = -rxE, 



(99) 



since Tq x ro = (remember that only for uniformly moving charge, E is directed along the present 
position ro). Both the fields go as and the Poynting vector vanishes over the surface at infinity; 
uniformly moving charge does not radiate. 

For an accelerated charge, it can be shown that (the deduction is not difficult but involves a 
lot of bookkeeping, and is done in any standard textbook; my favourite is Griffiths, which does a 
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brute force job) 

^ " 47r72r3(l - sin^ V)^/^ + 47rr3(l - J sin^ ^)3/2 't'' ' 

B = -r X E. (100) 
r 

The second term is the component that carries energy away to infinity (both E and B go as 1/r 
and hence the Poynting vector is nonzero even at infinity) ; this is the radiation field and exists only 
if a 7^ 0. Note that these fields are transverse: 

r.Erad = r.Brad = 0. (101) 

Therefore, 

S = E X B = - (E'^r - (r.E)E) = -E'^r. (102) 

7.1 Radiation from a Slow- moving Charge 

For a nonrelativistic motion, r tq, ?; — 0, so 

^rad = X (r X a) = £:[r(r.a) - a] (103) 

and 

where 9 is the angle between r and a. Thus, (i) there is no radiation along the direction of 
instantaneous acceleration; (ii) the radiation pattern is symmetric for 6 ^ tt — 6,Tr + 6 and 27r — 0, 
so shaped like a doughnut (an 8-shaped pattern) around the direction of a. 

To get the total power, we have to integrate S over the surface of the sphere with radius r: 



/ 



S.da = / sin2 ed(cos 6) ( d<P = ^.-.27r = —q^a^. (105) 

loTT^ J J 167r^ 3 DTT 



This is known as the Larmor formula, and is valid only for nonrelativistic motion. The power 
radiated does not depend on the sign of o; so whether the charge is accelerated or decelerated, we 
will get a radiation. 

Q. How should eq. (105) look in the SI and the Gaussian systems? [Hint: Don't just put back 
the factor of cq; make a dimensional analysis to see how many powers of c you need.] 



7.2 Relativistic Generalisation of Larmor's Formula 

How does eq. (105) change if the motion is relativistic? The deduction is due to Lienard, but we 
will follow a shorter path, using the covariance argument. Before that, let us just work out two 
derivatives: 

^7 d 2\-l 12 s s 

— = —(1 — V ) ' =7 v.a = 7 vacosa, 

at at 



dt 

where a is the angle between v and a, so that 

2 /j..\2 



/d(7v) 

\ dt J ydtj 



= 7^a^+7^i;^a^ cos^ a = j^a^{l+^'^v^)-^'^v'^a^ sin^ a = 7^ (^a^ - v^a^ siv? . 

(107) 
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Writing eq. (105) as 



P = 



dp dp 

GTrm^ V dt dt 



we immediately see that the relativistic generalisation is 

Qirm? \ dr dr 



(108) 



(109) 



where r = t/'j is the proper time, i.e., the time kept by the clock of the charge. Note that only 
dpn/dr, and not dp^/dt, is a Lorentz vector. Now 



dpfj, dpi^ 
dr dr 



dr. 



dE\' 



(110) 



and in the nonrelativistic hmit, t ^t, and the variation of E is negligible (since that is overwhelm- 
ingly controlled by the rest mass m), so we get back eq. (105). However, in the relativistic case, 
E = 7m, p = 7mv, and with eqs. (107) and (110), the expression for P reads 



P 



7 



dp\2 /dE^' 
Hit) ~\dt) 



(111) 



This is the relativistic generalisation of eq. (105). A comparison shows that there is a huge boost 
factor of 7^, so the emitted power is tremendously enhanced. Is this enhancement isotropic or is 

there a directional bias? This question we will address soon, but before that, let me just quote a 
formula for the angular distribution of radiated power, without deriving: 



dP |n X {(n — v) X a}p 

lin " 167r2 (1 - n.v)5 ' 

where n is the unit vector along r. An angular integration recovers eq. (111). 
Q. What should be the form of eq. (112) in the SI system? 



(112) 



7.3 Relativistic Motion: v || a 

Suppose that at a particular retarded time tr the velocity was instantaneously parallel to the 
acceleration. Then tq x a = r x a, and the electric field is given by eq. (103), except that the 
denominator will now have instead of r^. The same is true for B. But 



S 1 — V COS 

so we expect the Poynting vector to have the same form as in eq. (104) but an extra factor of 
(1 — vcosO)^ in the denominator. 

Now there is a catch. There will be a further correction of s/r on the energy loss, so that the 
denominator contains (1 — vcosO)^. Let us try to understand the origin of this correction, which 
is analogous to the well-known Doppler shift. 

The energy emitted by the electron in a time dt' is located in the volume between two spheres, 
one of radius r centred at x'2 and the other of radius r + dt' centred at x'l . Consider an infinitesimal 
volume element dv of this asymmetrical shell; suppose this subtends a solid angle dO, = dS/r'^ at 
x'2. Since dr = dt' — [r.\/r]dt', 

dv = dS dr = dS (l - dt' = ^dS dt'. (114) 
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Figure 3: Location of energy radiated by an electron as it moves from x'l to x^- 



Therefore the energy contained in this volume within the solid angle d£l is \{E'^ + B'^){s/r)dSdt' = 
E'^{s/r)dSdt'. 

Thus, the power radiated per unit solid angle, dP/dfl, is given by 

dP _ a? sin^ 6 

'dn ~ 167r2 {l-vcos9)^' ^ ' 

To get the total power, we have to integrate over Q. The integration over cp gives 2n, and the 
integration over 9 can be performed substituting x = cos 9 and using 

^ {l-x^)dx ^4 1 

-1 {l-vxf 3(l-t;2)3' ^-"^"^ 
so that 

Stt 3(1-^2)3 67r ' ^ ^ 

consistent with the Lienard formula, eq. (111). 

The power is not only enhanced, it is sharply peaked in the forward direction (see fig. 4). 
Though the power at precisely = is zero, most of it is concentrated within a narrow cone. The 
angle 9max where the radiated power is maximum can be obtained by differentiating dP/dO, with 
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12 



Figure 4: Plot of dP/dO. for w = (the 8-shaped pattern), v 
power has been scaled down by a factor of 10~^ for v = 0.98/ 



0.5 and v = 0.98. Note that the 



respect to 6 and setting the derivative to zero (of course, you have to check the second derivative 
to ensure that this is a maximum) . This gives 3f cos^ Omax + 2 cos Omc 



cos 



-1 



+ 1-1 



3v 



cos 



-1 



1 

8^ 



hv = 0, so 
1 

= 2^' 



(118) 



where the last fraction is at the limit v ^ 1. 



This is, very crudely speaking, what happens when a fast-moving electron decelerates (maybe 
within some material). The emitted radiation, called bremsstrahlung, is peaked in the forward 
direction. There is a frequency distribution of the spectrum; the maximum frequency of the emitted 



photon is obviously the change in the kinetic energy of the highest velocity electron, 



There is no lower limit on the frequency. There can be infinite number of zero-energy photons (the 
total energy is still zero!) emitted along the direction of v. 

Q. Show that at the limit u — > 1, Omax = I/27. (Hint: express v in terms of 7 and do a binomial 
expansion in terms of 1/7^.) 
Q. Find 

dmax for electrons with energy 1 GeV. 
Q. Assume both v and a to be along the z-direction: v 
the spherical polar coordinate: 



vk, a = ak. Let n be the unit vector in 



n = sin 6 cos 4>i + sin sin + cos Ok. 



(119) 



Get eq. (115) from eq. (112). 



7.4 Frequency Distribution: Bremsstrahlung for Slow Electrons 

Let us now go into a bit more detail of what has been said in the last paragraph of the preceding 
subsection. For simplicity, we will discuss the case v || a but w ^ 1, so that the retarded position 
and time are almost identical with the present position and time (s ~ r). 

Recall that in this limit, the electric field that contributes to radiation is given by eq. (103) and 
the Poynting vector by eq. (104). The amount of energy 8 flowing out per unit time, or the power 
radiated P, into a solid angle dVL, is given by 

d£ dP 9 ,^ ,99 g^a^sin^^ , „, 

|S|r2 = lE^rflV^ = ^ ^ . . (120) 



dtdn dQ ' ' ' 167r2 
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What is the spectral composition of the radiation? In other words, what is the amount of energy 
in the frequency band uj to cu + dto? For that, we first need to know how E^^d varies with the time 
t. Let us use the notation Et = \I^rad\{t)- We define the Fourier Transform (FT) of Et as 

1 foo 1 /■oo 

E^ = ^ Ete^^'dt, Et = ^ E^e-^^'dw. (121) 
y I'K J -oo yZ/KJ-oo 

We will also need Parseval's theorem, which says 

/•oo f-OO 

\Etfdt= \E^\'^duj = 2 lE^l'^dco, (122) 



where we use the last step to avoid dealing with negative frequencies, which does not have any 
meaning. 

From eq. (120), we have 

^ = 1 \Etfr'dt = 2r' J \E^fdw (123) 

so that 

df f 

2r'^\E^\^ ^ d£ = 2duj / \E^\^dS, (124) 



du dO. 

where we have used r'^dVL = dS. To go further you have to know the form of Et. But we can say 
something important even qualitatively. Recall that the FT of a narrow Gaussian (width~ a) is a 
broad Gaussian (widths l/cr) — this helps us to understand, for example, the uncertainty principle 
in quantum mechanics. So, if the velocity change occurs for a small time r (which defines, in some 
sense, the width of Et), the frequency spectrum will have a high-frequency cutoff at oJmax ~ l/r. 

Suppose the velocity change An of the charge takes place in a very short time interval At. If 
radiation takes place at toi Et during this time interval is proportional to Au/At, which we can 
write as AuS{t — to). This is nothing but making the velocity change instantaneous, and hence the 
time for the charge to interact with the retarding electric field zero. Thus, 



/oo 
iidt = Au. (125) 
-oo 



(Remember that all coordinates associated with the charge are actually retarded coordinates, but 
since the velocity is small, this distinction is irrelevant.) 

With such a simple ansatz, we have 



1 (126) 
V2^ 47rr ' ^ ' 



so that, from eq. (124), 



dS = 2 I lEj^dSdu 



Js 

= 2 ^(Au)^duj / d(j) / sin^^d^ 

2-K IbTT^ Jo Jo 



= ^duj{Au)\ (127) 

This means that the energy spectrum is uniform: equal amount of energy is emitted at every 
frequency interval, and there is no cutoff on lo. This is clearly unnatural, and occurs only because 
we have taken the interaction to be of zero duration. Actually, the energy of a bremsstarhlung 
photon cannot be greater than the kinetic energy of the charge, so 

TT^max = ^mu^. (128) 
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The number density of emitted photons in the frequency interval du is given by 

dN^ = T- = «^ Au — , 129 

where a = 1/137 is the fine structure constant. This tells you that as u goes to zero (photons 
becoming softer), the number density increases, and ultimately blows up for a; = 0. Thus, it is 
possible to shake off infinite number of zero-energy photons. This, obviously, does not affect the 

spectrum. 

Q. Show that in the SI system, the last equation in (124) will have a multiplicative factor of cqc 
on the right-hand side. 

Q. Show that in eq. (127), the right-hand side is to be multiplied by l/(eoc^). 

Q. Using cq. (128), draw the u versus £ plot. Redraw the same, replacing to by the wavelength A. 

Remember that du oc dX/X^. 

Q. Suppose that the velocity change is not instantaneous, but occurs for a time interval to — ^/'^ 
to to + during which the velocity changes uniformly and the total change is Atx. Show that 

l^^.l'oc^, (130) 

where ^ = uit/2. Plot the functional form of l-Et^jp. Do you have any idea why this looks like a 
diffraction pattern? 

7.5 Relativistic Motion: v ± a 

Suppose V and a are instantaneously perpendicular (for a circular motion, they are always perpen- 
dicular). Taking v = vk, a = ai, and eq. (119), we get 



nx{(n— v)xa} = a {—sin ^sin -|- cos ^(v — cos ^)}i -|- sin ^ sin ^ cos — sin ^ cos ^(i; — cos 

(131) 

Squaring, we get (after some trigonometric bookkeeping), using eq. (112), 



dP_ _ q^a^ [{l-vcos9y - (l-v^) sin^ 6> cos^ (p] 
'dn ~ 167r2 (1 -t; cos 6^)5 

An angular integration yields (do this) 



(132) 



^=^9'aV, (133) 

which you can get in a straightforward way from eq. (111). 

The radiation is sharply peaked in the forward direction along v, i.e., 6 = 0. Near 9 = 0, the 
0-dependent term is small, so the pattern is almost uniform in (p. For an electron in a circular 
orbit, the radiation sweeps like a beacon. This is known as the synchrotron radiation, after the 
machines where it was first observed. A charged particle revolving in a circular orbit loses energy, 
and as the rate of change of momentum dp/dt = ^yma, we have 

P = ^^-rH'^]\ (134) 



Btt V dt 

so that the energy loss is much more for an electron than a proton, for two reasons: m is smaller, 
and for the same energy, 7 is larger. This is precisely why the LEP is in all probability the last 
circular e^e~ collider, though protons will soon start running in the same machine with a much 
higher energy, but with significantly smaller synchrotron loss. (This machine is an engineering 
marvel: a perfectly circular pipe, deep underground, with a circumference of 27 km!) 
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The same phenomenon occurs in the sky. Pulsars are rotating neutron stars whose magnetic 

axis does not coincide with the rotational axis (the same is true for earth too) , and therefore relative 
to the fixed axis of rotation, the magnetic dipole vector is changing, and the star gives off magnetic 
dipole radiation. With a typical radius of 10 km, a rotational period of 10~^ s, and a surface 
magnetic field of 10^ Tesla, pulsars radiate away a huge amount of power in a highly directional 
beacon, and when that beacon intercepts the earth (it may not, in that case we won't be able to 
see it) we sec the pulsar. (This is a naive picture of a pulsar. On and near the surface of the 
pulsar, such radiations get damped because of the electromagnetic plasma present there, but the 
plasma in turn radiates, so we get some beacon-like pulsed radiation based on more or less the 
same principle.) 

Even for a circular motion, and more so for arbitrary motions, the frequency spectrum of the 

synchrotron radiation contains a large number of components. It can be shown (see Jackson, sec. 
14.4) that the highest frequency in the synchrotron spectrum, tOc, is given approximately by 

Uc ~ ^07^5 (135) 

where ujq is the frequency of rotation. Note the factor of 7^: in a typical 10 GeV machine, 
Jmax ~ 2 X 10"^, Wo ~ 3 X 10^ s^"*^, SO ujc ~ 2.4 X 10-*^^ s^^, corresponding to 16 keV X-rays. 

Q. Do the necessary vector algebra and trigonometry to get eq. (132) from eq. (112). 

Q. The last run of LEP had 105 GeV electrons. What was the value of 7maa;? What should be 

Imax for 7000 GeV protons? 

Q. Calculate the ratio of synchrotron losses (power radiated) for LEP-II (electrons with 105 GeV) 
and LHC (protons with 7000 GeV). Given, nip = 1836me. 



7.6 Thomson Scattering 



Suppose a plane monochromatic wave falls on a free particle of charge q and mass m. This will 
transfer some energy to the particle, so the particle will be accelerated and emit its own char- 
acteristic radiation. This radiation will be emitted in all directions (but will not be spherically 
symmetric). In the nonrelativistic limit (quantum-mechanically, where the energy of the photon is 
much less than the rest energy of the particle), however, the frequency of the scattered wave is the 
same as the incident wave. This process is known as Thomson scattering. 

In the nonrelativistic limit, eq. (112) reads 

The acceleration is provided by the incident plane wave. If its propagation vector be ko, and its 
polarisation vector ej„, the electric field can be written as 

B = einEoeM-ik^x'') (137) 

where k = (a;,ko). Remember that ej„ is a vector of unit magnitude and ej„.ko = 0, since the 
wave is transverse. For our subsequent analysis, we will take kg along the lab-fixed z-axis, so that 
without any loss of generality, two linearly independent choices of ei„ may be taken along the x 
and the y axes. 

The force equation for nonrelativistic motion is ma.{t) = qE, so 

a(t) = ein-Eo eM-ikf.x"). (138) 
m 

The scattered wave propagates along n (see figure), which makes an angle 9 with ko. If its polari- 
sation be e*„( (which can have two possible values, ei and 62), a little vector algebra shows that 



Since the polarisation vector is in general complex, we should complex conjugate the polarisation of the outgoing 
wave compared to the incoming wave. 
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(convince yourself of this!) 



so that 



|nx (nxa)p = |e*„t.a|' 



dP 



167r2' 



^out 



.a 



(139) 



(140) 



The electric field, and hence a, is a rapidly oscillating function of time, and we need an average 
over a complete cycle of oscillation. For nonrelativistic motion, a, though a function of time, does 
not change much over the oscillation period, and so the time-averaging is essentially computing 
\^in-^out\'^ and replacing the harmonic part by 1/2: 




Figure 5: Polarisation of the scattered wave. 
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The differential scattering cross-section, da/dft, is defined as the ratio of energy radiated per unit 
time (i.e., power) per unit solid angle to the incident energy flux per unit area per unit time. But 
the denominator is nothing but the time-averaged Foynting vector (the time-averaging brings a 
factor of 1/2), which is \Eo\'^/2, so 



da 

dn 



47rm 



l^m-^out I 



(142) 



The scattering geometry is shown in fig. 5. The polarisation vector ei lies in the plane containing 
n and ko, and 62 is perpendicular to it. In terms of unit vectors of the cartesian system i, j, and 
k, we may write 



ei = cos 6 cos 4)1 + cos 6 sin cf)] — sin ^k, 
62 = — sin (f)\ + cos (f)]. 



(143) 



In general we do not see the individual final-state polarisations, so we have to sum over all possible 
polarisation states. If the initial beam is polarised along the x-axis, the angular distribution. 
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summed over final-state polarisations, is (cos^^cos^^ + sin^^), while if it is polarised along the 
y-axis, the angular distribution is (cos^ ^sin^ cf) + cos^ <p). If the incident beam is unpolarised, we 
have to average over these contributions, so 



This is a procedure which is true for all scattering calculations, whether classical or quantum- 
mechanical: to get the scattering amplitude for unpolarised incident beam, average over all initial- 
state polarisations and if you do not look for a specific polarisation in the final state, sum over all 
final-state polarisations. 

Eq. (144) is the Thomson formula for scattering of electromagnetic radiation by a charged 
particle. It is valid only if the energy of the incident radiation is much less than the rest energy of 
the particle, e.g., scattering of waves upto X-ray (~ keV) by electrons or upto 7-ray (~ MeV) by 
protons. The distribution is symmetrical around 9 = tt/2 (and does not have any dependence), 
and by integrating, we get the Thomson cross-section: 



For an electron (or a proton), /A-k is the fine-structure constant a (in the natural system), which 
is approximately 1/137. The ratio a/m, for an electron, is approximately 2.82 x lO""-*^^ m. This is 
called the classical electron radius. Of course, electron is a point particle, but a classical distribution 
of charge, which total the electronic charge, must have a radius of this order if its electrostatic self- 
energy is equal to the electron mass. Note that the spin of the particle does not play any role in 
the scattering 

Q. In the natural system of units, a = 1/137, and electron mass m = 0.511 MeV. Also, 197 MeV-fm 
ft! 1 (energy and length have opposite dimension, right?). Prom this, show that a^/m ft 2.82 fm 
= 2.82 X lO-i^m. 

Q. What should be the formula for the classical electron radius in the SI system? 

Q. Calculate the Thomson cross-section for electrons in barn (1 barn = 10~^^ cm~^). 

Q. Suppose the electron is a uniformly charged sphere of radius tq, the classical electron radius. 

If its angular momentum is h/2, calculate the maximum linear velocity of a point on its surface. 

What's wrong with this value? 

Q. Suppose we look at only those scattered waves which are polarised along ei. Show how the differ- 
ential Thomson cross-section, da/dVl, will vary with 6 for a fixed value of ^ (take = 0, 7r/4, 7r/2, vr) 
and with ^ for a fixed value of (take = 0, 7r/4, 7r/2, tt). Repeat the exercise with waves polarised 
along 62. 

Q. If we sum over the final-state polarisations but start with polarised incident wave, show how 
da/dCl will look like as functions of 9 and (f) for (i) beam polarised along x-axis; (ii) beam polarised 
along y-axis. Take the same values of 6 and as given in the previous problem. 

7.7 Modifications to the Thomson formula: Compton, Klein-Nishina 

When the energy of the incident radiation is comparable to the rest energy of the particle, mod- 
ifications to eq. (144) occur. The calculation is quantum-mechanical in the sense that it uses the 
concept of photons, but the result can be translated into the language of classical electrodynamics. 

^^This is expected since spin docs not have a role in classical physics. Quantum-mechanically, it is rather subtle to 
show that the same cross-section is obtained in the nonrelativistic limit for spin-0 and spin-1/2 particles. However, 
there is a difference: spin-1/2 particles have a nonzero magnetic moment which interacts with the wave and contributes 
to the scattering matrix element, see the next subsection. 




(144) 




(145) 
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Considering the scattering of a photon of energy E = huj hj a, particle of rest mass m and applying 
the two-body relativistic kinematics, we get the Compton formula: 

^' ' (146) 



E 1 + A(i_cos0)' 



where E' is the energy of the scattered photon (or the wave, if you prefer), and 6 is the scattering 
angle. This modifies the unpolarised differential cross-section: 



Compton Thomson \E 

The modification factor comes entirely from the phase-space, and da/dVL drops at large 9 from its 
Thomson value, the drop depending on the ratio Ejm. To be precise, eqs. (146) and (147) are true 
only for scattering from spin-0 particles. To get the total cross-section, one must remember that 
E' /E is also a function of ^. At limiting cases, we get 

f E ^ 2>m , „^ , ^ 

— = 1-2— {E<^m), -— m<E . 148 
or rn ' AE ^ 

For spin-1/2 particles like electrons, not only the charge but also the magnetic moment is responsible 
for scattering. This modifies the matrix-element squared part of eq. (147): 

le- e* J2 >\e- e* + ^ fl49) 

where k and k' are the wave numbers of the incident and the scattered waves respectively. This is 
known as the Klein-Nishina formula. 

Q. Plot the differential scattering cross-section as a function of 9 for (i) Thomson; (ii) Compton 
with E/m = 0.5, and (iii) Klein-Nishina with E/m = 0.5. 

Q. Show that for E <^ m, the Compton cross-section is indeed given by a = aril — 2E/m). 



7.8 Scattering from Bound Electrons 

Consider an electron bound to an atom, with a characteristic frequency loq. This is, of course, 
a semiclassical description, where we assume the electron to be bound in a harmonic oscillator 
potential, so that without any external force, the equation of motion is 

mx -I- mujQ-x. = . (150) 

If we apply a harmonic electric field E = einEo exp{—iLot) on this electron, the equation of motion 
will be 

mx + tuujqX = q'E = qeinEo exp{—iujt) . (151) 
The solution is (which should be familiar to you if you have studied optical dispersion) 

yi = ein— i^T- — 9 ) exp{-iujt) , (152) 

m \loq — w^/ 

which can be checked by direct substitution in eq. (151). The acceleration, in turn, is given by 

a = X = Bin— I n ^ ^ ) exp(-za;t) . (153) 
m Kuq — J 

For free electrons, ujq 0, which is the Thomson scattering that wc have just studied. For strongly 
bound electrons, lvq^ uj, and da/dO., which is proportional to a^, is scaled by a factor of {lo/loq)^. 
This is known as Rayleigh scattering; we can write 

(^Rayleigh = ( ) (^Thomson ■ (154) 
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Thus, the cross-section is large for higher frequencies. This explains why the rising and setting suns 
are red; the rays have to cross a thicker layer of atmosphere and bluish components are scattered 
away. This also explains why the sky is blue; if we are not in the line of the sun, there is more 
chance to catch the scattered blue components than the red ones, which are less scattered. The 
red range of the visible spectrum is at 4100 A while the violet end is at 6500 A; thus, we expect a 
(41/65)^ 16% scattering of the violet end compared to the red end. 

What happens when loq = ojl This is when the electrons vibrate with their resonance frequency, 

the displacement is maximum, and so is the absorption. The refractive index, instead of increasing 
slowly with w, drops sharply; this is known as the region of anomalous dispersion. In the next 
section, we will see how the pole at w = is physically removed. 

7.9 Cherenkov Radiation 

We have said that a uniformly moving charge does not radiate. There is one exception. Note the 
factor (1 — vcosO)^ in the expressions of E and B. Since v < 1, this factor is always positive. 
But suppose a charged particle is moving in a medium where its velocity exceeds that of light 
in that medium. The denominator will now read (1 — micosO)^, where n is the refractive index 
of the medium, and can be zero \i 6 = cos~^(l/nt;). So there will be, theoretically, an infinite 
field intensity in the direction specified by 9: a cone about the velocity vector. Experimentally, 
one gets radiation inside this cone; this is known as Cherenkov radiation. This is nothing but an 
electromagnetic shock wave where the surface of the cone is the shock front. Outside the cone 
the field intensities are zero; on the surface they arc infinity and inside the cone they are finite but 
nonzero. So if a highly energetic particle enters a medium, we may get a light signal at an angle 9. 
This effect was first observed by Frank and Tamm, and the theoretical explanation was provided 
by Cherenkov. 

The Cherenkov effect is used extensively in high-energy particle detectors. Note that if you 
can successfully track the emitted light (with photon detectors like photomultiplier tubes), you can 
accurately determine the velocity, and hence the energy, of the particle. For cosmic-ray and neutrino 
experiments (like SuperK), where the direction of the incoming particle is extremely important to 
know, Cherenkov detectors are essential. 

The frequency distribution of the emitted radiation is continuous but the intensity increases 
with increasing frequency. That is why most of the Cherenkov radiation is in the ultraviolet; even 
the visible radiation appears blue. 

Q. Two electrons, with energies 5 MeV and 10 MeV respectively, enter water, whose n = 1.33. 
Calculate the Cherenkov angles. 

8 Radiation Reaction 

An accelerated charge radiates, so it loses energy. Thus, we expect that under a given force, its 
acceleration will be smaller than a particle with same mass but no charge. This is as if the emitted 
radiation is exerting a back force on the charge. This is known as radiation reaction. 

So far, we have never included this effect in our calculations. For example, in the calculation 
of bremsstrahlung or synchrotron radiation, we have started from the expression of E^a^, and 
calculated the Poynting vector, and then integrated the flux over a surface at infinity. We have 
never considered the fact that velocity and acceleration may change due to the reactive effects of 
radiation. 

^^An acoustic shock wave is characterised by abrupt change of pressure and temperature; they always increase 
inside the shock, and the boundary is known as the shock front. 
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Hopefully, we can justify ourselves; the experiments match with calculations. Thus, the radia- 
tion reaction should be small in such cases. It is better that this is so; because, as we will see soon, 
this is one case which classical electrodynamics is just not competent enough to deal with, and so 
is one of the major motivations for a quantum theory of electromagnetic fields. 



8.1 When is the Radiation Reaction Important? 

There are, in general, two types of motion that are physically interesting. First, the steady accel- 
eration, where the particle accelerates uniformly from t = to t = T, and its velocity changes from 
zero to aT. Second, the periodic motion, where the particle vibrates with a characteristic frequency 
ujQ and an amplitude d. Let us consider both these cases one by one. 

Case 1: Steady acceleration — The radiated power is given by the Larmor formula: P = 
q^a^/Gir. Note that this is a non-relativistic formula; everything that wc derive now is non- 
relativistic (a relativistic generalisation will be considered later). Thus, the total energy loss in 
time T is of the order of 

2 2 

Erad=^T. (155) 
DTT 

The typical energy of the system, Eq, is of the order of 

Eo ~ mv'^ = m{aTf . (156) 

If Erad (not to be confused with the radiation electric field) is to be comparable with £'o, we should 
have 

= ma^T^ ^ T = ^ . (157) 

67r Dvrm 

This value of T is known as the characteristic time and will henceforth be denoted by r. Thus, 
radiation reaction is significant if T r and negligible if T S> r. 

Let us try to have an estimate of r for electrons. We use a = /Air = 1/137, mass of the 
electron m = 0.511 MeV, and 197 MeV-fm = 1. This gives r = 1.876 x 10~^^ cm = 6.25 x 10"^^ s. 
This is indeed a very small time; typically, light takes this much time to cross an atomic nucleus. 

Case 2: Periodic motion — With an amplitude d and a characteristic frequency ljq, the typical 
energy Eq ^ miOQd^, while the acceleration a ~ cogd and the time period T ~ I/coq. The energy 
radiated in a full cycle is (g^a^/67r)T, which should be compared with the typical energy Eq. 
Putting the relevant quantities, we get the condition for radiation reaction to be significant over a 
full cycle: 

cjQT ~ 1 . (158) 

Thus, we can neglect the radiation reaction if r ^ "^o^^- This is true, for example, even in the most 
exotic cases like fast-rotating pulsars. However, the radiation effects can build up over time, and 
the pulsar should fractionally slow down over a million of years. 

Thus, we conclude that the radiation reaction is important only if (i) the motion changes 
appreciably over a time scale r, or (ii) over a length scale cr (in our system, c = 1, so both should 
be of the order of r). This justifies the neglect of this effect in our earlier calculations. 

What is the energy scale associated with the time scale r? For electron, this is just 1/r = 105 
MeV. An electron with such an energy is highly relativistic, while all the derivations that we are 
going to perform are non-relativistic in nature. Thus, the radiation reaction is always a small 
perturbation, if it is at all significant. 
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8.2 The Abraham-Lorentz Formula 



The rate at which the charged particle loses energy is given by the Larmor formula: 

2 2 

F.ad.v = -^, (159) 

where the minus sign indicates that the energy is being lost. However, this equation is incorrect. 
The energy is carried away to infinity by the radiation fields; but there are other components of 
E X B, viz., those that go as and r~^, which carry energy but do not take it far away from 
the configuration. When the charge undergoes a periodic motion, energy is pumped back and forth 
into these components, and they exert a back force as well All is not lost, however, since if 
we take a time interval at the end of which the system returns to its previous state, the net force 
exerted by the and components vanish. Thus, we would be safe if we write 

Frad-^ dt = -^ a.a dt, (160) 

Jtl DTT Jtl 

with the proviso that v(ii) = \{t2), a(ii) = a(t2)- One can integrate the right-hand side of eq. 
(160) by parts: 



whose first term vanishes by the periodic boundary conditions, so that 

rt2 
iti 

Thus, the simplest possible form of Fj.ad can be written as 

F™d = |^a. (163) 

This is known as the Abraham-Lorentz formula. However, note that this is not a true derivation; 
we do not know, for example, the component of F^.^^ perpendicular to v. 

An exact derivation is a bit cumbersome. One has to assume a charge distribution; this, in 
itself, is a big objection, as for all we know the electron is a point charge. Without going into the 
derivation, let us just note down the main features. 



Jtl 



Frad ~ ^a 
DTT 



.V dt = 0. (162) 



• The charge must be distributed over a small but finite region, as the field due to a point charge 
blows up on the charge (E — *■ oo as r — 0). In his original derivation, Lorentz assumed a 
spherical charge distribution. One can have other distributions; e.g., Griffiths considers a 
dumbbell-shaped configuration. 

• The next step is to break the charge distribution into several fragments. For an accelerated 
charge, the force exerted by fragment 1 on fragment 2 is not the same (and opposite) to that 
of fragment 2 exerted on fragment 1. This is a breakdown of Newton's third law, but can be 
justified from the fact that the retarded position of 1 with respect to 2 is not the same as the 
retarded position of 2 with respect to 1. The imbalances, summed up, give the net radiation 
reaction, or the force of the charge on itself. 

^^This is why Griffiths calls Radiation reaction a misnomer; according to him it should be called Field reaction. 
However, the former term sticks. 
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• The first term of the expression for the self-force is a term that is proportional to acceleration, 
and thus just adds up to the 'bare' mass of the charge; this is the mass that wc will observe 
in an experiment. There is a problem; in the limit when the distribution is shrunk to a 
point, this term blows up. Thus, if we consider the electron as a point particle, the mass 
would have been infinite! Clearly, this is a very fundamental problem, related to the exact 
nature of electron. We need a quantum theory to address this issue, but let me mention that 
the infinities haunt the quantum theories as well. Of course, there we use different tricks to 
tackle this problem. And we could do that very successfully, as is evident from the success of 
quantum electrodynamics. 

• To circumvent this problem, Lorentz assumed the bare mass of the elctron to be zero; the 
mass that we observe is the result of this self-force. It can be shown that the radius needed 
for the charge distribution in this case is exactly the classical electron radius. However, we 
have seen that such a classical picture is untenable; the velocity of a point on the 'equator' 
of electron is far above the velocity of light! 

• The second term in the expression for the self-force is the Abraham-Lorentz term, proportional 
to a. Rest of the terms drop out in the limit of a point charge. So we can work with the 
Abraham-Lorentz formula in the classical limit. But, all in all, this is not a success story. 



8.3 Problems with the Abraham-Lorentz Formula 

Now for the paradoxes. Suppose the acceleration is constant: a = 0. The charge radiates, but the 
radiation reaction is zero! 

The opposite extreme, is a charge that is momentarily moving with constant velocity: a = 
but a 7^ 0. The charge does not radiate, but there is a self-force acting on it! 

Consider the motion of a charge without any external force. The equation of motion is 

Fr-ad = g;^a = ma , (164) 

whose solution is 

a(f) = a{t = 0) exp(t/T) , (165) 

where r is the characteristic time for the charge. Unless a.{t = 0) is zero, this is a runaway solution: 
the acceleration increases spontaneously with time! 

One can, however, try to circumvent the runaway solution. Suppose there is an external force, 
which can be written as a function of time, so that one can write 



Let us try a solution 
so that 

and 

Integrating this, we get 



mv - mrv = Fext{t) . (166) 
v(i) = exp(i/T)u(i) , (167) 

v(i) = - exp(t/r)u(t) + exp(t/T)u(i) , (168) 
r 

mu{t) = --eM-t/r)'Pext{t) . (169) 
r 

1 /■* 

mu(t) = / exp{-t'/T)Fextit') dt' , (170) 

r Jc 
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where C is the integration constant. This can be written as 



T Jc 



(171) 



Thus, acceleration at any time t depends not only on the instantaneous value of the external force 
but sort of its weighted time-average. In particular, the acceleration at t will depend on the value 
of the force at some earlier time t'\ If t — < 0, the preacceleration is significant only for a time of 
the order of r, but it is there: even before the force is applied, the particle starts accelerating! This 
violates causality, but at the microscopic scale: r is of the order of 10~^^ s. Still, this is repugnant 
from a theoretical point of view. One can circumvent either the runaway solution or the acausal 
preacceleration, but not both. 

Is there such a preacceleration? In other words, can an electron feel a force before it is applied? 
This is difficult to answer from an experimental point of view: such preacceleration is significant 
only for a very small time which is not accessible even to the present-day experiments, and thus 
macroscopic causality is not violated. Theoretically, this would need, in contrast to the retarded 
potentials discussed earlier, an advanced potential, a signal propagating in the backward light-cone. 
Definitely, we will be more comfortable with a theory which does not involve this aesthetically 
repelling component. 

And these problems do not go away even in a relativistic classical theory. 



8.4 Relativistic Motion: Dirac Formula 



A Lorentz-covariant generalisation of the equation of motion would be 
and the covariant generalisation of {q'^/6Tr)k can be written as 



(172) 



Si 

Gtt 



(173) 



where S is some scalar. Now, the energy radiated per unit time, which is the zero-th component 
of Pi^u^ divided by 7, is equal to the work done per unit time, F^ad-v, so that energy conservation 
implies 

(174) 



But u'^U/j, = 1, so this gives 



S = —u. 



d^uf __d_i' dur\ du^ dui" 
ds \ ^ ds J ds ds 



The first term, which is ^d{unU^)/ds, is zero, as u^U/j, = 1, so 

„2 



(175) 



d'^u'^ „ dui, dvy 
|- 1 

ds' ds ds 



(176) 



This is known as the Dirac formula. Note that the last term is proportional to the expression for 
power in the relativistic generalisation of Larmor's formula: 



P = - 



dun du^ 
Gtt ds ds 



(177) 
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8.5 Radiation Reaction on a Charged Oscillator 



Consider a charged particle, typically an electron, bound in a harmonic oscillator potential with a 
characteristic frequency lvq. Since the charge is accelerating all the time, one may wish to include 
the radiation reaction too for the equation of motion. Just for simplicity, let us consider the 
one-dimensional case. The equation of motion is 

mx + mcoQX — Frad = , (178) 

which can be written as 

x + ujqX-tv = 0, (179) 

as Fj-ad = TTiTV. Without the last term, the solution would have been x{t) = xq exp(ia;oi); the 
particle vibrates with a frequency ivq. Let the solution be of the form 

x{t) = xoe^p{-at) (180) 

where a can be a complex number. One obtains, from eq. (179), 

TQ^ + a^+w^ = 0, (181) 

whose solution, for r = 0, is the well-known a = ziziujQ. The cubic equation can be exactly solved; 
let us just note that there is one real solution, the so-called runaway solution, where the acceleration 
goes on increasing, and hence the solution is unphysical. The remaining two solutions are complex, 
conjugate to each other, since both r and ujq are real. In the limit ujqt <C 1, we just quote the 
result from Jackson, which is true upto 0{u!qT^): 

a = ^r±(wo + Au;), T = lo^t , Alo = —cu^t^ . (182) 

Note that the solution has an exponentially falling term exp(— |rt) and an oscillatory term. This is 
reminiscent of a damping force that is proportional to x; the result will be similar whenever we have 
odd number of time derivatives of x. The frequency has been shifted by Alu ((this is known as the 
level shift), which is indeed very small, unobservable to experiments, since Auj/uq ~ {ooot)^ ^ 1- 
The decay constant T causes the energy of the oscillator, which goes as x^, to fall with an exponential 
pattern exp(— Ft). Classically, this should be the width of a spectral line. 

However, this classical treatment is not applicable to electronic states in an atom since the 
spectral lines are quantum mechanical in origin. One can perform a more sophisticated calculation 
in the ambit of quantum field theory. The qualitative results, the existence of a nonzero level shift 
and a radiative broadening, are there, but their magnitudes differ from the classical counterparts. 
In particular, the level shift can be much larger, even larger than the line width (classically, the 
former is always smaller than the latter by a factor of u)qt). One famous example is the shift of 
hydrogen 251/2 level; without any radiative corrections, this should be degenerate with the 2Pi/2 
level, as we expect from the Dirac theory of electrons. Experimentally, the 25i/2 level is about 
1000 MHz lower than 2Pi/2j this is known as the Lamb shift after Willis Lamb, who discovered it 
(with his student Retherford), and is one of the cornerstones of quantum electrodynamics. The 
Lamb shift is entirely due to the radiative effects. 

Q. Consider the equation of motion of a charged oscillator in the presence of a damping, so that 
the equation of motion is 

x + ^x + ulx = . (183) 

Assuming a solution x{t) = a;oexp(— at), show that a is indeed complex. Find the decay width 
and the level shift if 7 ^ wq- 
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8.6 Scattering and Absorption of Radiation by a Charged Oscillator 



Next consider a charge in a harmonic oscihator potential with characteristic frequency ujq being 
driven by a driving force of frequency ui. Suppose this is an electric field; for one-dimensional case, 
I can write 

X + loqX — TV = — \Eo\ exp(—iijt) . (184) 
m 

The system oscillates with frequency 00, so d^x/dt^ = — w^i, and 

X + ^x + ulx = —\Eq\ , (185) 
m 

where 7 = uP't. The solution is similar to that of eq. (152): 

X = — \Eo\ \ ) exp(-za;t) . (186) 

m \uJq — uj'^ — iLO^ J 

Since 7 ^ cj, we can neglect the radiation reaction both in the Thomson limit [loq ^ oj, almost 
free electron) and in the Rayleigh limit {uj 3> wq, strongly bound electron). When w « (Jq, the 
energy does not have a pole any longer; it has been flattened out to a Gaussian shape by the oj^ 
term (note that in this limit, 7 = uj'^t = ujqT = F of eq. (182)) This is the range where the 
dispersion is anomalous and the absorption is strong, just because the oscillator vibrates with a 
large amplitude, and can take up lots of energy from the source. 



^'A similar term occurs in the quantum theory when we consider the arnphtude for a particle to move from some 
point A to some other point B. The amplitude has a pole when the particle is real, i.e., it satisfies the energy- 
momentum relationship p'^p^ = w? . This pole is removed by a similar term if the particle is unstable and can decay. 
This is just a reflection of the uncertainty principle: for a particle which is not stable, the energy, or the rest mass, 
cannot be measured with absolute precision. The form of the amplitude with such a pole-removing term is known as 
the Breit-Wigner form. 
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